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Crystal structure, electronic structure, and magnetic behavior of the spin-| quantum magnet 
Nai.5VOP04Fo.5 are reported. The disorder of Na atoms leads to a sequence of structural phase 
transitions revealed by synchrotron x-ray powder diffraction and electron diffraction. The high- 
temperature second-order a fi transition at 500 K is of the order-disorder type, whereas the 
low-temperature /3 -f-)- 7 + 7' transition around 250 K is of the first order and leads to a phase 
separation toward the polymorphs with long-range (7) and short-range (7') order of Na. Despite 
the complex structural changes, the magnetic behavior of Nai.5VOP04Fo.5 probed by magnetic 
susceptibility, heat capacity, and electron spin resonance measurements is well described by the 
regular frustrated square lattice model of the high-temperature a-polymorph. The averaged nearest- 
neighbor and next-nearest-neighbor couplings are J\ ~ —3.7 K and J2 ~ 6.6 K, respectively. Nuclear 
magnetic resonance further reveals the long-range ordering below Tjv = 2.6 K in low magnetic 
fields. Although the experimental data are consistent with the simplified square-lattice description, 
band structure calculations suggest that the ordering of Na atoms introduces a large number of 
inequivalent exchange couplings that split the square lattice into plaquettes. Additionally, the 
direct connection between the vanadium polyhedra induces an unusually strong interlayer coupling 
having effect on the transition entropy and the transition anomaly in the specific heat. Peculiar 
features of the low-temperature crystal structure and the relation to isostructural materials suggest 
Nai.5VOP04Fo.5 as a parent compound for the experimental study of tetramerized square lattices 
as well as frustrated square lattices with different values of spin. 

PACS numbers: 75.50.Ee, 75.30. Et, 75.10.Jm, 61.66.Fn 



I. INTRODUCTION 

Among a broad family of low-dimensional and frus- 
trated magnets, systems based on the frustrated square 
lattice (FSL) model enjoy special attention from the- 
ory and experiment. The model entails competing ex- 
change couplings along sides (Ji) and diagonals (J2) 
of the square and is mostly studied for purely Heisen- 
berg Hamiltonian and spin-^. Extensive theoretical 
worl5pJ^ convincingly established three ordered ground 
states that emerge for different values of the frustra- 
tion ratio a — J2I J\- The Neel and columnar antifcrro- 
magnetic (AFM) states are separated by a critical spin- 
liquid region around the quantum critical point at a = ^ 
(Refs.[3]-[H]). At a = — ^, the columnar AFM phase should 
border the ferromagnetic (FM) phase. However, the pre- 
cise nature of this boundary remains controversial. While 
an earlier report proposed a nematic phase separating 
the regions of FM and columnar AFM ground states 
Richter et a/.* demonstrated the single abrupt transition 
without any intermediate phases around a = — \. 

Experimentally, the FSL-type magnetic behavior has 
been observed in V"'"'* compounds with [VOXO4] layers 
comprising VO5 square pyramids and non-magnetic XO4 
tetrahedrgPttS (the only known exception is PbVOa with 
magnetic layers formed by VO5 pyramids exclusiveljl^^ . 



If X is a main group (p) cation, the leading coupling is 
AFM J2, whereas J\ is usually weaker and can be ei- 
ther FM_(2C= P) or AFM (X = Si, Ge), see also Ta- 
ble [yn p ^ l ^'^ l^ Recent efforts in crystal growtff^ and 
experimental investigatiorPsHlo] gycj^ c ompou nds es- 
tablished the columnar AFM ground statd^IltSSI and re- 
vealed simple trends in thermodynamic properties j ^^ l ^^^ l ^^ -^ 
in line with theoretical predictions.^ A thorough compar- 
ison between experiment and theory, however, spots cer- 
tain discrepancies. For example, the sublattice magne- 
tization is gradually reduced from 0.6 /is in Li2VOSi04 
(a ~ 10) to 0.5 [IB in Pb2VO(P04 )2 (a ^ -2) and 
0.4 /is in SrZnVO(P04)2 [a ~ -l).EiE2l By contrast, 
theory predicts a nearly constant sublattice magnetiza- 
tion of 0.6 [iB in this range of a (Ref.[H|. The discrepancy 
may arise from the spatial anisotropy of the FSL due to 
the low crystallographic symmetry of Pb2VO(P04)2 and 
SrZnVO(P04)2. In fact, none of the reported FSL-type 
phosphates, ^A'VO(P04)2 {AA! = Pb2, PbZn, SrZn, 
BaZn, BaCd) with FM Ji, are tetragonal, thereby the 
effects of the spatial anisotropy are expected.l^ 

To facilitate the experimental verification of theo- 
retical results for the FSL model, tetragonal systems 
with a perfect square lattice of magnetic atoms are 
required. Motivated by this challenge, we explored 
the structure and properties of Nai.5VOP04Fo.5. This 
compound was prepared in 2002 by a hydrothermal 
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TABLE I. Lattice parameters, space groups, and refinement residuals for diflerent polymorphs of 
Nai.5VOP04Fo.5. 



Polymorph 

a 
P 
7 
7' 



Temperature 
(K) 
560 
298 
150 
150 



(A) 

6.39563(1) 
9.03051(2) 
12.76716(2) 
6.37996(4) 



(A) 

10.65908(2) 
10.62002(3) 
10.57370(4) 
10.5910(1) 



Space group 

lA/mmm 
P4:2/mnm 
P42/mbc 
lA/mmm 



Ri/Rp 

0.031/0.109 
0.035/0.103 
0.033/0.073 
0.022/0.073 



method. Massa et al^ found a tetragonal crystal struc- 
ture with the three-dimensional (3D) VOPO4F0.5 frame- 
work formed by the FSL-type VOPO4 layers (Fig. [T]). 
Sauvage et alW^ claimed to prepare the same compound 
by a high-temperature annealing in air, and reported a 
very similar structure refinement. Since Na atoms in 
Nai.5VOP04Fo.5 and related syst erns a re readily dein- 
tercalated and exchanged with Li,l^32ll the frustration 
ratio a could be tuned, making Nai.5VOP04Fo.5 an ap- 
pealing FSL system. To explore this possibility, we per- 
formed a comprehensive study of the parent compound 
Nai.5VOP04Fo.5. 

The outline of the paper is as follows. We list experi- 
mental and computational procedures in Sec.llTl and pro- 
ceed to details of the crystal structure in Sec. |III[ Further 
on, we perform thermodynamic and magnetic resonance 
measurements (Sec. IV I and band structure calculations 
(Sec. |V]) to evaluate the low-temperature magnetic be- 
havior and the microscopic magnetic model. The com- 
parison of Nai.5VOP04Fo.5 to known FSL compounds 
and other structural analogs is given in Sec. |VI| followed 
by a summary and outlook. 



II. METHODS 

Powder samples of Nai.5VOP04Fo.5 were prepared by 
a solid-state reaction of Na4P207, VO2, and NaF in 
an evacuated and sealed silica tube at 700 °C for 24 
hours. The stoichiometric mixtures of reactants were 
pelletized, placed into corundum crucibles, and cov- 
ered with lids to avoid the reaction between silica and 
NaF. Na4P2 07 was obtained by the decomposition of 
Na2HP04 in air at 400 °C. The bluish-green powders of 
Nai.5VOP04Fo.5 were single-phase, as confirmed by lab- 
oratory x-ray diffraction (XRD) measured with Huber 
G670 Guinier camera (CuKqi radiation, 29 = 3 — 100° 
angle range, image-plate detector). 

High-resolution XRD data for structure refinement 
were collected in the 150 — 560 K temperature range 
at the ID31 beamline of European Synchrotron Radia- 
tion Facility (ESRF) with a constant wavelength of about 
0.4 A. The signal was measured by eight scintillation de- 
tectors, each preceded by a Si (111) analyzer crystal, in 
the angle range 29 — 1 — 40 deg. The powder sample 
was contained in a thin-walled borosilicate glass capil- 



lary with an external diameter of 0.5 mm. The sample 
was cooled below room temperature in a He-flow cryo- 
stat and heated above room temperature with a hot-air 
blower. To achieve good statistics and to avoid the effects 
of the preferred orientation, the capillary was spun dur- 
ing the experiment. The JANA2006 program was used 
for the structure refinement.!^ Symmetry changes at the 
structural phase transitions were analyzed with ISODIS- 
PLACE program.lSfi- 

The samples for an electron diffraction (ED) study 
were prepared by crushing the powder in ethanol and 
depositing the suspension on a holey carbon grid. ED 
patterns were taken at room temperature (RT) and at 
100 K using a Philips CM20 microscope equipped with a 
Gatan cooling holder. 

The differential scanning calorimetry (DSC) measure- 
ment was performed with a Perkin Elmer DSC 8500 in- 
strument in the temperature range 120 — 800 K in ar- 
gon atmosphere with a heating/cooling rate of 10 K/min. 
The powder sample of Nai.5VOP04Fo.5 was placed into 
a corundum crucible. 

The magnetic susceptibility of Nai.5VOP04Fo.5 was 
measured with an MPMS SQUID magnetometer in the 
temperature range 2 — 380 K in applied fields up to 5 T. 
Heat capacity measurements on a pressed pellet were per- 
formed by relaxation technique with a Quantum Design 
PPMS instrument in the temperature range 0.45 — 200 K 
and in fields up to 11 T. The data below 1.8 K were 
collected with a ^He insert. 

The electron spin resonance (ESR) measurements were 
performed at X-band frequency {v = 9.36 GHz) on 
a Bruker ELEXSYS E500-CW spectrometer, equipped 
with a continuous He-gas flow cryostat (Oxford Instru- 
ments) operating in the temperature range 4.2 — 300 K. 
The polycrystalline powder sample was fixed in a quartz 
tube with paraffin and mounted in the center of the mi- 
crowave cavity. The field derivative of the microwave- 
absorption signal was detected as a function of the static 
magnetic field due to the lock-in technique with held 
modulation. Resonance absorption occurs when the in- 
cident microwave energy matches the energy of magnetic 
dipolar transitions between the electronic Zeeman levels. 

The nuclear magnetic resonance (NMR) measurements 
were carried out using pulsed NMR techniques on ^"'^P 
(nuclear spin I = ^ and gyromagnetic ratio 7Ar/27r = 
17.237 MHz/T) nuclei in the temperature range 1.5 — 
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TABLE II. Atomic positions, isotropic atomic displace- 
ment parameters (f/iso, in 10^^ A^), and occupancy fac- 
tors (/) of a-Nai.5VOP04Fo.5 at 560 K (upper lines) and 
7'-Nai.5VOP04Fo.5 at 150 K (bottom lines). The / values 
are fractions of total occupancy for a given position. 



FIG. 1. (Color online) Crystal structure of 

a-Nai.5VOP04Fo.5: the layers in the ab plane (left) 
and the projection along c (right). The layers have an 
ideal FSL geometry with the couplings Ji and J2 (left) and 
connect along c via the F atoms (right). Disordered Na 
atoms occupy the voids of the resulting framework (right). 



300 K. We did the measurements at radio frequencies of 
79 MHz and 8.62 MHz, which correspond to applied fields 
of about 4.583 T and 0.5 T, respectively. The spectra 
were obtained either by sweeping the field or by doing the 
Fourier transform. The NMR shift K = (iJ^of - H) /H 
was determined by measuring the resonance field of the 
sample (H) with respect to a nonmagnetic reference 
H3PO4 (resonance field i?rcf)- The '^^P spin-lattice relax- 
ation rate 1/Ti was measured by the conventional single 
saturation pulse method. 

Scalar-relativistic band-structure calculations for 
Nai.5VOP04Fo.5 were performed within the framework 
of density functional theory (DFT) using the basis set 
of local orbitals (FPL09. 01-35 code).'2ij We applied the 
local density approximation (LDA) with the exchange- 
correlation potential by Perdew and Wang'^^ and a k 
mesh of 216 points in the first Brillouin zone (40 points 
in the irreducible wedge). Exchange couplings were 
evaluated by mapping V Sd bands onto a multi-orbital 
Hubbard model (see Sec. |v] for further details). 



III. CRYSTAL STRUCTURE 
A. High-temperature a-polymorph 

According to RT studies of Refs. US! and HZl 
Nai.5VOP04Fo.5 has a body-centered tetragonal unit 
cell with flsub ^ 6.38 A and c ~ 10.62 A. Our RT XRD 
pattern is largely consistent with this unit cell, although 
10 weak reflections remained unindexed. These reflec- 
tions can be assigned to the v^Osub x V^Osub x c su- 
percell and evidence the superstructure formation at RT. 
The superstructure reflections in Nai.5VOP04Fo.5 disap- 
peared upon heating. Above 500 K, the patterns could be 
fully indexed on the Osub x ftsub x c body-centered tetrag- 
onal unit cell that we further refer as a-modification, in 
contrast to /3-modification at RT (Table |l]). 

The structure refinement •^■^ for the a-phase (Table Hi]) 



Atom Position 


X 


y 




^ ISO 


f 
./ 


Nal 


8h 


0.2651(5) 0.2651(5) 





2.9(2) 


0.52(1) 






0.2627(8) 0.2627(8) 







0.49(2) 
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161 


0.394(3) 
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11(1) 






0.336(3) 
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0.22(6) 
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1 















0.22(6) 


1 



conforms to the atomic positions given in Ref. [26l We 
find a 3D framework of corner-sharing VO5F octahedra 
and PO4 tetrahedra (Fig. [T]) . Vanadium atoms form one 
short bond to the 01 atom (along the c direction), four 
longer bonds to oxygens in the ab plane (the 02 position), 
and one long bond to the fluorine atom (also along c), see 
Table Hn) The octahedra are linked via the PO4 tetrahe- 
dra in the ab plane, whereas the fluorine atoms connect 
the octahedra into pairs along c. Na atoms occupy the 
voids of the resulting framework. 

At this point, the question concerning the arrangement 
of O and F atoms may arise. Since atomic numbers of 
these elements differ by unity, it is hard to distinguish 
between O and F using XRD, especially in a powder ex- 
periment. The single- crystal refinementP^ proposed the 
complete ordering of the O and F atoms. The ordered 
arrangement of these atoms is also supported by empiri- 
cal arguments. The 01 position corresponds to the short 
V-0 bond which is typical for oxy gen and uncommon 
for fluorine, e.g., in oxyfluorides.'^^ Further on, the 02 
atoms in the ab plane are parts of rigid PO4 tetrahe- 
dra and can not mix with fluorine. Therefore, fluorine 
is left to its 2a position with two long bonds to vana- 
dium. Note also that our ESR and NMR experiments 
(Sec. 



IV) suggest unique (or very similar) positions for 
vanadium and phosphorous, respectively, thus indicating 
the ordered arrangement of O and F atoms below RT. 

While the O and F atoms in Nai.5VOP04Fo.5 form an 
ordered framework, the Na atoms are disordered. Our 
refinement for the a-modification identifies two Na po- 
sitions, 8h (Nal) and 16^ (Na2). After constraining the 
sum of occupancies to the Nai 5VOPO4F0.5 composition 
and refining atomic displacement parameters (ADPs) to- 
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TABLE III. Interatomic distances (in A) in a-, and 7-polymorphs of Nai.5VOP04Fo.5 at 560 K, 
298 K, and 150 K, respectively. 





a 




P 




7 


V-Ol 


1.626(3) 


V-Ol 


1.636(2) 


V-Ol 


1.629(3) 


v-02 


4 X 2.011(1) 


V-02 


2 X 2.001(5) 


V-02 


2.032(9) 






V-03 


1.991(5) 


V-03 


1.986(9) 






V-04 


2.023(5) 


V-04 


1.992(9) 










V-05 


2.01(1) 


P-02 


4 X 1.530(1) 


Pl-02 


4 X 1.543(6) 


Pl-02 


4 X 1.549(8) 






P2-03 


2 X 1.540(7) 


P2-03 


4 X 1.542(8) 
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2 X 1.509(7) 
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2 X 1.502(8) 
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gether with the occupancies, we found that the 8h posi- 
tion is approximately half-filled, whereas the 161 position 
is filled by | (i.e., twice less Na atoms than in 8h).^^ This 
is in line with Ref. 26 that reported the occupancies of 
I and ^ for 8ft. and 16^, respectively (our experiment 
is done at 560 K compared to 300 K in Ref. (26', hence 
a redistribution of Na atoms is possible). By contrast, 
Ref. 127] assigns Na2 to a 8j position with a negligible oc- 
cupancy of jq, i.e., 0.5 atoms per unit cell in 8j compared 
to approximately 2 atoms in 161 in our refinement. 



B. Room-temperature /3-polymorph 



The superstructure formation in the RT /?-polymorph 
is confirmed by ED (Fig. [2|. While the intense re- 
flections are assigned to the body-centered tetragonal 
flsub X flsub X c unit cell of the a-polymorph, weaker su- 
perstructure reflections are clearly visible in the [001], 
[010], and [101] patterns. XRD and ED suggest a prim- 
itive tetragonal v^Csub x V^f^suh x c unit cell, with 
the reflection condition hOl, h + I — 2n that identifies 
the P4:2/mnm space group. This refiection condition 
is present in the [010] pattern (Fig. [2]). The emergence 
of the forbidden reflections, such as 100 and 010 in the 
[001] pattern, is caused by the multiple diffraction. Fi- 
nally, the [110] pattern reveals sharper reflection condi- 
tions h + k = 2n and h + I ^ 2n corresponding to the 
body-centered unit cell of the a-polymorph. The lack of 
superstructure reflections, such as 110, in the [110] ED 
pattern is due to their zero structure factors, in agree- 
ment with the complete absence of these reflections in 
the XRD pattern. 

The superstructure formation at room temperature is 
related to the partial ordering of the Na atoms in t he 
/3-modification. The structure refinement (Table IV l'^'^ 
showed that the VOPO4F0.5 framework remains intact 
(Table [ni| ), whereas Na atoms are partially ordered in 
two 8i positions. The Nal position is completely filled, 
and the Na2 position is exactly half-occupied under the 
constraint of the Nai.5VOP04Fo.5 composition. 



FIG. 2. ED patterns of /3-Nai.5VOP04Fo.5 at RT. The pat- 
terns are indexed on a tetragonal v^Osub x \/2asub x c su- 
percell. Bright reflections correspond to the flaub x a^uh x c 
subcell, whereas weak reflections indicate the supercell. 



TABLE IV. Atomic positions and isotropic atomic displace- 
ment parameters (;7iso, in 10"^ A^) of /?-Nai.5VOP04Fo.5 at 
298 K. The Na2 position is half-filled. 



Atom 


Position 


X 


y 


z 




Nal 


8i 


0.5144(5) 


0.2314(7) 





2.8(1) 


Na2 


8i 


0.785(1) 


0.028(1) 





6.0(4) 


V 


8j 


0.2484(1) 


0.2484(1) 


0.80118(6) 


0.60(1) 


PI 


id 





1 

2 


1 
4 


0.6(1) 


P2 


4e 








0.2558(3) 


0.2(1) 


01 


8j 


0.2499(4) 


0.2499(4) 


0.3529(2) 


1.3(1) 


02 


16fc 


0.0957(6) 


0.4031(6) 


0.1623(6) 


0.7(2) 


03 


8j 


0.0949(6) 


0.0949(6) 


0.1662(7) 


0.5(2) 


04 


8j 


0.4037(6) 


0.4037(6) 


0.1617(8) 


1.0(3) 


F 


4/ 


0.2450(4) 


0.2450(4) 





0.8(1) 
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FIG. 3. (Color online) Temperature dependence of 
the subcell volume (top) and c/osub ratio (bottom) for 
Nai.5VOP04Fo.5. Error bars are smaller than symbols. Lines 
are guide for the eyes, with arrowheads showing changes upon 
heating/cooling in the region of the /3 7 + 7' transition.'^ 
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TABLE V. Atomic positions and isotropic atomic displace- 
ment parameters ((7iso, in 10~^ A'^) of 7-Nai.5VOP04Fo.5 at 
150 K. 



Atom 


Position 


X 


y 


z 




Nal 


8h 


0.3715(4) 


0.6291(4) 





0.69(5 


Na2 


8h 


0.1186(4) 


0.5951(4) 





0.69(5 


Na3 


8h 


0.3486(4) 


0.3764(3) 





0.69(5 


V 


16i 


—0.0021(1) 


0.2511(1) 


0.1991(1) 


0.35(3 


PI 


Ah 
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4 


0.29(3 


P2 


4d 





1 
2 


\ 

4 


0.29(3 


P3 


8<7 


0.2490(2) 


0.7490(2) 


1 

4 


0.29(3 


01 


16i 


-0.0144(3) 


0.2440(4) 


0.3522(3) 


0.28(4 


02 


l&i 


0.0957(7) 


0.0079(4) 


0.3400(6) 


0.28(4 


03 


l&i 


0.5965(7) 


-0.0044(4) 


0.3379(6) 


0.28(4 


04 


16i 


0.2583(3) 


0.3476(7) 


0.3261(6) 


0.28(4 


05 


16i 


0.7567(4) 


0.3440(8) 


0.3436(6) 


0.28(4 
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8h 


-0.0045(4) 


0.2463(4) 





0.28(4 
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FIG. 5. (Color online) XRD pattern of Nai.5VOP04Fo.5 
collected at 150 K (black dots) and the two-phase refinement 
(solid line). Upper and lower ticks denote the reflections of 
the 7'- and 7-phases. The inset shows the unit cells of a- 
(solid), /3- (dotted), and 7- (dashed) polymorphs. 



FIG. 4. (Color online) DSC data showing the first-order 
/3 -f-)- 7 -I- 7' transition at 265 K (heating, upper panel) and 
235 K (cooling, bottom panel). 



The structure refinements for a- and /3- 
Nai.5VOP04Fo.5 suggest that the phase transition 
at 500 K is of the second order (order-disorder type). 
Indeed, the temperature evolution of lattice parameters 
(Fig. |3]) shows a smooth change in the cell volume around 
500 K. The structural change also conforms to sym- 
metry requirements for a second-order transition: the 
P4:2/mnm space group can be derived from M/mmm 
using the irreducible representation. Similar to 
second-order transitions in (CuCl)LaNb207 (Ref. l37| . 
we did not observe the a O /3 transformation by DSC, 
presumably, due to the small change in the entropy. 



C. Low-temperature phase separation 

Below RT, one can expect a further ordering of the 
Na atoms. Indeed, we found another structural trans- 
formation around 250 K. This transition is of the first- 
order type, as shown by: i) the temperature hysteresis 
(about 235 K on cooling and about 265 K on heating. 
Fig. [4]); ii) the coexistence of the low-temperature and 
high-temperature phases in a certain temperature range; 
iii) the abrupt change in the cell volume (Fig. [s]). The 
transition is also revealed by an abrupt change in the 
ESR Hnewidth (see Sec.|lVC|. 

Upon cooling the sample, new reflections appeared be- 
low 235 K, and the reflections of the /3-phase fully dis- 
appeared at 210 K. Below 210 K, the pattern could be 
indexed on a 2asub x 2asub x c tetragonal unit cell that we 
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FIG. 6. ED patterns of 7-Nai.5VOP04Fo.5 at 100 K. The 
patterns are indexed on a tetragonal 2asub x 2asub x c unit cell. 
Brighter reflections with even h, k, and I indices arise from 
the subcell, whereas weaker spots evidence the formation of 
the superstructure. 
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FIG. 7. ED patterns of 7'-Nai.5VOP04Fo.5 at 100 K. 
The patterns are indexed on a body-centered tetragonal 
Ssub X flsub X c unit cell. Diffuse lines in the [110] pattern 
evidence short-range order. 



further refer as 7-inodification. However, a closer exam- 
ination showed a complex shape of the (hkl) reflections 
with non-zero /. In particular, the (002) reflection at 
29 = 4.34° could not be fitted as a single peak (Fig. [5]). 
This points to the coexistence of two phases with slightly 
different lattice parameters. 

The low-temperature ED study confirmed the presence 
of two different phases below RT. One of the phases is 
the ordered 7-polymorph that shows sharp superstruc- 
ture reflections of the 2asub x 2asub x c unit cell (Fig. |6]) . 
The crystallites of the second phase are less ordered and 
reveal diffuse satellites or diffuse intensity lines instead of 
the superstructure reflections (Fig.]?]). The 7-polymorph 
has a primitive tetragonal unit cell. The reflection con- 
ditions hhl, I = 2n ([110] pattern) and hOl, h = 2n ([001] 
and [010] patterns) were confirmed by XRD, and resulted 



in the PA2/rabc space group. Similar to the [110] pattern 
of /3-Nai.5VOP04Fo.5 (Fig. [2|, the [010] pattern of the 
7-polymorph contains subcell reflections only. 

The second low-temperature phase, further referred as 
7'-polymorph, reveals sharp subcell reflections as well as 
diffuse satellites on the [001] pattern and diffuse lines 
along c* on the [110] pattern (Fig.[7|. The intensity of the 
diffuse lines is modulated in a way that the maxima co- 
incide with the superstructure reflections of the 7-phase. 
This suggests that the two low-temperatures polymorphs 
develop a similar supercell, with the long-range order in 
7 and the short-range order in 7'. While the ordered 
structure of 7-Nai.5VOP04Fo.5 can be determined from 
XRD (see below), the nature of the short-range order in 
7'-Nai.5VOP04Fo.5 is difficult to establish because the 
compound is unstable under the electron beam. 

In XRD, the diffuse intensity produced by the 
7'-polymorph is smeared. The remaining subcell reflec- 
tions are formally equivalent to that of the a-phase. The 
150 K pattern could be refined as a two-phase mixture 
of the 7-phase with the 2asub x 2asub x c unit cell and 
the 7'-phase with the Osub x agnh x c unit cell. The re- 
fined composition of the mixture yields about 40 % of 
the 7'-phase. This ratio is nearly constant within the 
temperature range under investigation. 

The structure refinement of 7-Nai.5VOP04Fo.5 (Ta- 
ble |v|'^'^ identifies three Na positions which are fully oc- 
cupied, as evidenced by their low ADPs.l^ The ordering 
of Na has moderate effect on the framework: the V-0 
and P-0 distances are nearly unchanged compared to 
the a- and /3-phases (Table III I. The refined structure of 
7'-Nai.5VOP04Fo.5 closely resembles the a-polymorph 
(Table [ll]). Based on the XRD and ED data, we sug- 
gest that Nai.5VOP04Fo.5 develops a two-phase mixture 
at the first-order transition around 250 K. Specific heat, 

(Sec. 



IV C and Ref. 



DSC, and ESR measurements 
do not reveal any further structural changes below this 
temperature. 

At low temperatures, the 7- and 7'-phases have the 
same subcell volume but a different c/flgub ratio (Fig. |3|. 
Based on our synchrotron XRD data as well as mag- 
netization, NMR, and ESR measurements (Sec. IV I, 
we conjecture that both phases contain V^**, and their 
composition is, therefore, similar (within the available 
resolution). The origin of the phase separation in 
Nai.5VOP04Fo.5 is presently unclear and should be dis- 
closed in future studies. For a further discussion on Na 
ordering in Nai.5VOP04Fo.5 and related compounds, we 
refer the reader to Sec. IVI Al 



IV. MAGNETIC BEHAVIOR 

A. Magnetization 

The temperature dependence of the magnetic suscep- 
tibility (x) of Nai 5VOPO4F0.5 is typical for a low- 
dimensional (and possibly frustrated) antiferromagnet. 
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see Fig. [s] At low fields (/xoff < 1 T), the maximum at 
T^g^x ~ 5.7 K manifests the onset of AFM short-range 
correlations followed by a kink at Tjv — 2.6 K due to the 
long-range ordering (see also the specific heat and NMR 
spin-lattice relaxation rate data in Figs. 10 and 16 re- 



spectively). A stronger magnetic field has a pronounced 
effect on x'- the maximum shifts to lower temperatures, 
whereas slightly increases (see also Fig. 11). Above 



50 K, l/x(T) shows a nearly linear behavior indicating 
the Curie- Weiss paramagnetic regime. No signatures of 
the /3 O 7 + 7' transition around 250 K could be ob- 
served. 

Above 50 K, the fit of the data with the Curie- Weiss 
law: 



X = Xo + 



c 



T + e 



(1) 



yields temperature-independent susceptibility xo = 
— 1.5(1) X 10"'* emu/mol (core diamagnetism and 
Van Vleck paramagnetism), the Curie constant C = 
0.366(2) emu K/mol, and the Curie-Weiss temper- 
ature 9 = 3.3(1) K. The effective magnetic mo- 
ment of 1.71(1) /is conforms to the expected value 
of g^B\/S{S + 1) — 1.69 for where we used 

and the powder-averaged g — 1.95 from ESR 



(see Sec. [rVC |. 

Despite the 3D nature of the crystal structure, a low- 
dimensional magnetic behavior should b e expected. Pre- 
vious studies of compound^^HSHini identify the typ- 
ical dxy orbital ground state driven by the short V-01 
bond along the z direction. Since the d^y orbital docs not 
overlap with the orbitals of the axial 01 and F atoms, 
exchange couplings along the c direction are weak. In 
the case of Nai.5VOP04Fo.5, this implies the formation 
of VOPO4 magnetic layers, similar to AA'YO(F04)2 
FSL-ty pe vana dium phosphates (see the left panel of 
Fig. nj) | ii | i4 | 2 4] 'pjji-ee different phosphorus positions in 
the 7-Nai.5VOP04Fo.5 structure lead to a number of 
inequivalent exchange couplings (Sec. [v| and Fig. 18). 
However, the thermodynamic properties in the high- 
temperature region (T > | Ji|, J2) should follow the pre- 
dictions for the ideal FSL model, and allow one to eval- 
uate Ji and J2, the averaged nearest-neighbor (NN) and 
next-nearest-neighbor (NNN) couplings on the square 
lattice.ISI 

Based on the above arguments, we fit the mag- 
netic susceptibility of Nai.5VOP04Fo.5 with the high- 
temperature series expansion (HTSE) for the FSL 
model!^ 



X = Xo 



knT 



Ji 



(2) 



where xo is similar to Eq. Q, Na is Avogadro's num- 
ber, g is the g-factor, and Cmn with m, n < 9 are 
coefficients of the HTSE. Owing to the internal sym- 
metry of the model, one typically finds two equivalent 
fits with different Ji and J2 (see Refs. [TTJ [TSl [Ml and 
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FIG. 8. (Color online) Temperature dependence of the 
magnetic susceptibility (x) and the HTSE fit with Eq. ([2]): 
Tmin = 12 K, Ji = -3.7(1) K, J2 = 6.6(1) K. The inset shows 
the field dependence of x low temperatures. The kink at 
2.6 — 2.8 K denotes the onset of the long-range order, see also 
Fig.[ll] 
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FIG. 9. (Color online) Details of the magnetic susceptibil- 
ity fit: averaged exchange couplings Ji and J2, g- value, and 
temperature-independent contribution xo evaluated for the 
fitting ranges with different Tmin- 



[25|) . For example, fitting the data above 12 K, we ar- 
rive at FM Ji = -3.7(1) K and AFM J2 = 6.6(1) K 
{g = 1.973(2), Xo = -1-5_(1) x lO""* emu/mol) or AFM 
Ji = 5.5(1) K and FM J2 = -2.1(1) K (5 = 1.981(2), 
Xo — —1.6(1) X lO"'^ emu/mol). The correct solution 
can be determined from the specific heat data^the 
saturation field,'231or band structure calculations .'^^'251 
the case of Nai.5VOP04Fo.5, band structure calculations 
provide a solid evidence for leading AFM couplings be- 
tween next-nearest neighbors (see Sec. Iv]). The analysis 
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of the saturation field also favors the solution with FM 
J I and AFM J2 (see below). 

A more subtle problem of the HTSE fit lies in the 
choice of the fitting range. The convergence of the 
HTSE depends on the frustration ratio a, hence the lower 
limit of the fitting range (Tmin) is not universal. Fit- 
ting the data with different Tmin, we find a stable so- 
lution for Tmin = 12 - 16 K (Fig. [9|. The resulting 
Xo = -1-5(1) X 10^4 emu/mol and g = 1.973(2) are 
in excellent agreement with the Curie- Weiss fit {xo = 
-1.5(1) X 10"'' emu/mol, g = 1.97(1)). At lower T^in, 
the high-temperature part of the HTSE fit deviates from 
the Curie- Weiss fit, thereby less accurate estimates of Ji 
and J2 are obtained. Using Ji ~ —3.7 K and J2 ~ 6.6 K, 
we arrive at a = J2/J1 — —1.8 that coincides with the 
frustration regimes of Pb2VO(P04)2, BaZnVO(P04)2, 
and PbZnVO(P04)2, see Table [VH] as well as Refs. |MI 
and [23 

The magn etizati on curve of Nai.5VOP04Fo.5 (right 



panel of Fig. 11 is typical for a two-dimensional (2D) 
antiferromagnet. The linear increase in the magnetiza- 
tion [M) at low fields is followed by a slight curvature 
and the saturation at /io-ffs — 15.4 T. To test our es- 
timates of Ji and J2 against iJg, we use the expression 
from Ref. SUfor the columnar AFM state (J2 > -0.5Ji) 
on the regular square lattice: 



/io-ff, = 2(Ji -t- 2J2)kB/{gfiB) 



(3) 



Assuming g — 1.95 from ESR (Sec. |IV C I , one finds 
^qHs ~ 14.5 T which is about 1 T lower than the ex- 
perimental value. The solution with AFM Ji implies the 
Neel ordering (antiparallel spins on nearest neighbors), 
leads to noH^ ~ 16.8 T, and overestimates Hs, which in 
this case is proportional to 4Ji instead of 2( Ji -I- 2 J2) in 
Eq. ([3]) ^ While the saturation field itself does not allow 
us to choose the correct solution, DFT results (Sec.Jvj) 
put forward the sizable AFM J2 and, therefore, FM Ji - 
AFM J2 regime with an additional interlayer coupling. 
This interlayer coupling contributes to the saturation 
field, and explains the underestimate of in the purely 
2D model described by Eq. By contrast, the AFM Ji 
- FM J2 regime overestimates Hg even in the 2D model, 
thereby any additional couplings will only exaggerate the 
discrepancy between the exchange integrals and satura- 
tion field. This way, we choose the FM Ji - AFM J2 
solution, which is confirmed b y the sim ilarity to other 
FSL-type vanadium phosphates.'^^'SlIIIl 



B. Heat capacity 

The measured specific heat {Cp}^- of Nai.5VOP04Fo.5 
resembles that of BaCdVO(P04)2 and other FSL 
compounds .'^^'21111! xhe sharp anomaly at Tjv — 2.6 K 
is followed by a maximum around 4 K and an increase in 
Cp above 8 K (Fig. 
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FIG. 10. (Color online) Zero-field specific heat (Cp) of 
Nai.5VOP04Fo.5 with a transition anomaly at Tjv = 2.6 K 
and the maximum of the magnetic contribution around 4 K. 
The inset shows the low-temperature data plotted as Cp/T 
against . Solid line indicates the idealized behavior. 



8 K are dominated by the phonon contribution. Unfortu- 
nately, we were unable to fit the high-temperature part 
with Debye functions and extract the magnetic contribu- 
tion. The fitting problem may be related to the coexis- 
tence of 7- and 7'-polymorphs having different phonon 
spectra. 

The specific heat maximum around 4.0 K suggests 
weaker exchange couplings in Nai.5VOP04Fo.5 com- 
pared to Pb2VO(P04)2 and PbZnVO(P04)2 that reveal 
J2 = 9 — 10 K and the maxim a of Cmag at a higher 
temperature of 4.5 - 5.0 K.li^By contrast, Li2VOSi04 
with a s imilar J2 = 6.3 K has the maximum of Cmag 
at 3.6 Since the phonon contribution at 4 K is 

small, the measured value of C„ 



the magnetic contri 



10 1. The maximum is a signature of 



jution Cmag, whereas the data above 



^ is a reasonable esti- 
mate for the maximum of Cmag- This way, we find 
Cmag/^ — 0.43 that compares well to 0.44, as expected 
for the unfrustrated^ or moderately frustratecP' square 
lattice, and expe rimentally observed in a number of FSL 
compoundsl ^^ l ^^ l 

The low-temperature evolution of Cp resembles clas- 
sical antiferromagnets. However, the data do not fol- 
low a simple behavior (see the inset of Fig. 10 1. 
The origin of this deviation is presently unclear. Note 
that other FSL-type compounds also show complex 
low-temperature features in the specific heat. The 
data for Pb2VO(P04)2 and Li2VOGe04 can be de- 
scribed by a combination of cubic and linear terms.f^^ In 
Nai.5VOP04Fo.5, the linear term is vanishingly small, 
but the cubic term itself poorly fits the data below 2 K. 

Heat capacity measurements in magnetic field es- 
tablish the temperature-vs-field phase diagram of 
Nai.5VOP04Fo.5 (Fig. 11). With increasing magnetic 
field, the specific heat maximum is suppressed, whereas 
the transition anomaly is growing and Tjv increases up 
to 2.8 K at 5 T. Above 5 T, both the anomaly and 
T/v are gradually suppressed. The overall shape of the 
phase boundary is typical for 2D antiferromagnets. ^^ ' ^^ ' ^ 
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FIG. 11. (Color online) Left: specific heat of Nai.5VOP04Fo.5 measured in different applied fields. Right: the field-vs- 
temperature phase diagram, based on the heat capacity and NMR measurements, as well as the magnetization isotherm at 
1.4 K (Ref.[ll|). 



The initial increase in Tjv is due to the field-induced 
anisotropy that weakens quantum fluctuations. Upon 
further increase in the field, the tendency towards the 
fully polarized state competes with AFM correlations and 
impedes the AFM long-range ordering. 

C. Electron spin resonance 

Typical ESR spectra of Nai 5VOPO4F0.5 are shown 
in Fig. |12[ The compound reveals a single exchange- 
narrowed resonance line which is perfectly described by 
a Lorentzian curve in the whole temperature range under 
investigation. The temperature dependence of the corre- 
sponding fit parameters, i.e., the resonance field -ffrcs, 
intensity /esr,; and half-width at half-maximum A_ff, 
is illustrated in Fig. |13| At room temperature, the g- 
factor, derived from the resonance field H-^cs — 3422 Oe 
via hv = gji^Hj-es is found to be g = 1.95, which is a 
typical value for (Sd^) ions coordinated by oxygen 
atoms in a square pyramid or an octahedron. On ap- 
proaching Tjv, the value decreases only slightly. The 
integrated ESR intensity nicely matches the static sus- 
ceptibility X obtained from SQUID measurements, and 
exhibits a maximum around 7 K. This indicates that all 
vanadium spins are observed by ESR. 

The linewidth is AiJ ~ 110 Oe at room temperature, 
and exhibits two anomalies upon cooling. One anomaly 
is close to Tjs[, probably due to short-range order effects, 
while the other one matches the structural /3 7 + 7' 
transition at Tgtr — 250 K. Similar to the structural 
phase transition in the quasi one-dimensional magnet 
CuSba 06,1221 the change in the ESR linewidth for T < T^tr 
is well described in terms of a thermally activated behav- 
ior: 

Ai/(T) = Ai/o + A.exp(-^) (4) 

with the BCS-like mean-field gap A(r) = 1.74Ao • (1 - 
r/rstr)"-^ We obtain Aq = 560(13) K, T^tr = 257(1) K, 



1 


1 ' 

4 K 

1 


i 1 i 1 i 

^^^^^^^^^^^^^^^^^ 


' 1 ' 

V = 9.36 GHz ^ 
1 


1 ' 

300 K 

1 



I 1 I 1 I 1 I 

3.0 3.3 3.6 3.9 

H (kOe) 



FIG. 12. (Color onhne) ESR spectra of Nai.5VOP04Fo.5 
measured at X-band frequency for selected temperatures in 
the paramagnetic regime. Solid lines indicate fits with the 
field derivative of a Lorentzian. 



A ~ 13(1) Oe, and the residual linewidth AiJo = 
93(1) Oe. It is interesting to consider the relation be- 
tween the energy gap Ag and the transition tempera- 
ture Tstj., which yields 2A/rstr = 4.36. This is enhanced 
compared to the value of 3.52 predicted by mean-field 



10 



3.48 



2 3.45 



3.42 



^0.30 
E 

E 

-2- 0.15 



0.00 



03 

i, 0.10 

X 

< 



0.08 



g = 1 .94 

,/ 



Na VPO F 

1.5 5 0.5 



9 = 1.95 

V'^(3d^ ) 



+ 



' \ ' 

o ESR 

• SQUID (H = 5 kOe) 




9.36 GHz 



100 



200 



300 



T(K) 



FIG. 13. (Color online) Temperature dependence of the ESR 
parameters: resonance field -ffros (upper frame), ESR intensity 
compared to the magnetic susceptibility (middle frame), and 
linewidth AH (lower frame) in Na1.5VOPO4F0.5- The solid 
red line indicates the fit using a mean-field BCS-like ansatz 
[Eq. (|4f]. 



theory for a second-order phase transition. Indeed, the 
transition at Tstr is of the first order, as shown by the 
temperature evolution of unit cell parameters (Fig. [s]). 
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FIG. 14. (Color online) Fourier-transformed ^^P NMR 
spectra at different temperatures T (T > Tjv) for 
Nai.5VOP04Fo.5. 
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FIG. 15. (Color online) Temperature dependence of the ''^P 
NMR line shift K measured at 79 MHz. The inset shows K 
plotted against the magnetic susceptibility x with tempera- 
ture as an implicit parameter. The solid line represents the 
linear fit with Eq. 



D. Nuclear magnetic resonance 



For an / = 2 nucleus, one expects a single spectral 
li)2 (j39 | 47 | 4 8i that is indeed observed in the NMR spectra 
of Nai.5VOP04Fo.5 in the 1.5 — 300 K temperature range 
(Fig. 14 1. Although our structural study (Sec. Ill) iden- 



tifies two inequivalent phosphorous positions at 300 K 
(/3-phase) and at least three inequivalent phosphorous 
positions in the 7-phase below 200 — 250 K, this will not 
necessarily lead to separate lines in the spectra. The 
line shift for an individual P site is determined by its lo- 
cal environment and the hyperfine couplings to electronic 
spins on V"*"*. Since the inequivalent P sites have very 
similar coordination and similar connections to the four 



surrounding vanadium atoms (Table III and Fig. [18| , we 
believe that the line splitting is too small to be observed 
in the present experiment. While this does not allow to 
refine our model of the low-temperature crystal structure 
and phase separation, the single spectral line facilitates 



the study of the magnetic properties by NMR. 

Since our measurements are done on randomly- 
oriented polycrystalline samples, the asymmetric shape 
of the spectra corresponds to a powder pattern due 
to an asymmetric hyperfine coupling constant and an 
anisotropic susceptibility. The line position was found 
to shift with temperature, similar to the magnetic sus- 
ceptibility (Fig. [Tsj ) . To establish the relation between K 
and x-i we use the expressionj2^ 



(5) 



where Kq is the temperature-independent chemical 
shift, is the hyperfine coupling constant, and 

Xspin/NA is the intrinsic (spin) magnetic susceptibility of 
Nai.5VOP04Fo.5 in units of /xs/Oe per electronic spin 
(i.e., per formula unit). In order to calculate ^hfi we 
plotted K vs. X with T as an implicit parameter (the 



11 



100 r 



10 r 




v^::^.ooo-oo-ooooo-o-o-o-oocQ^ 



_A— 8.62 MHz (//|^= 0.5 T) 
— O— 79 MHz 4.583 T) 



10 100 
Temperature (K) 

FIG. 16. (Color online) Temperature dependence of the spin- 
lattice relaxation rate 1 /Ti measured at two different frequen- 
cies of 79 and 8.62 MHz. 
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SrZnVO(P04)2 from ^ip NMR (Ref.|2ni), VOM0O4 from 
95Mo NMR (Ref. and Cu(HC02)2 ■ 4D2O from 
NMR (Ref. in]). The decrease in 1/Ti can be explained 
by the cancellation of antiferromagnetic spin fluctuations 
at the probed nuclei. 

The sharp peak in l/Ti identifies the onset of the long- 
range magnetic ordering. The peak position depends on 
the NMR frequency and reveals the field dependence of 
T/v. In Fig. [11] we plot the two transition temperatures 
derived from our NMR experiments, and find perfect 
agreement with the phase boundary obtained from the 
heat capacity measurements. 

Below Tjv, the ^^P line measured at 79 MHz broad- 
ens abruptly. In order to check whether any extra fea- 
tures could be resolved, we remeasured the spectra at a 
lower frequency of 8.62 MHz (Fig. [T7|. Two symmetrical 
shoulder-like features develop on both sides of the cen- 
tral line and move systematically away from each other 
with decrease in temperature. The symmetric position 
and the systematic evolution of the shoulders in the low- 
frequency low-temperature spectra is an indication of a 
commensurate magnetic ordering. Our low-temperature 
spectral shape can be compared to thejLi-NMR spectra 
of Li2VOSi04 reported by Melzi et alW^In the columnar 
AFM phase, they observed a central line and two sym- 
metrical satellites from the single-crystal measurement. 
In our case, such satellites appear as broad shoulders due 
to the random distribution of the internal field for the 
measurements on polycrystalline samples. Nevertheless, 
the overall spectral shape is consistent with ^Li-NMR re- 
sults on Li2VOSi04 and hence points to the columnar 
AFM ordering in Nai.5VOP04Fo.5. 



FIG. 17. (Color online) Temperature-dependent ^^P NMR 
spectra measured at 8.62 MHz. 



MICROSCOPIC MAGNETIC MODEL 



inset of Fig. The hnear fit with Eq. ^ over the 

whole temperature range yields ^hf = (212 ± 20) Oe/^s 
that is comparable to hyperfine couplings found in other 
FSL-type vanadium phosphates.!^ 

For an / = I nucleus, the recovery of the longitudinal 
magnetization is expected to follow a single-exponential 
behavior. In Nai.5VOP04Fo.5, the recovery of the nu- 
clear magnetization after a saturation pulse was indeed 
fitted well by the exponential function 



1 - 



M{t) 
Mo 



(6) 



where M{t) is the nuclear magnetization at a time t after 
the saturation pulse, and Mq is the equilibrium magneti- 
zation. The temperature dependence of 1 /Ti is depicted 
in Fig. 16 Above 10 K, 1/Ti is essentially temperature- 
independent, which is typical for paramagnetic moments 
fiuctuating fast and at random.'SSl Below 10 K, l/Ti de- 
creases slowly, shows a peak around 2.8 K, and further 
decreases at low temperatures. A similar behavior has 
been observed in Pb2VO(P04)2 from ^^P NMR (Ref.,23|). 



Despite the tetragonal crystallographic symme- 
try, the complex low-temperature superstructure of 
7-Nai.5VOP04Fo.5 gives rise to a number of inequivalent 
superexchange pathways and the overall distortion of the 
FSL. Using the atomic positions for 7-Nai.5VOP04Fo.5 
(Table |v| , we arrive at six different couplings in the ab 
plane: Ji, J(, and Jf between nearest neighbors as well 
as J2, ^^'^ ^2 between next-nearest neighbors, see 
Fig. 18 The spin lattice entails regular Ji — Jj' pla- 
quettes alternated with less regular plaquette-like units 
having one type of the NNN couplings (either J2 or J2) 
yet different NN couplings. 

While it is hardly possible to evaluate all the six afore- 
mentioned parameters experimentally, DFT calculations 
are capable of providing microscopic insight into indi- 
vidual superexchange pathways and establishing realistic 
spin models for complex crystal structures. "^^ In the 
following, we present DFT results for the ordered struc- 
ture of 7-Nai.5VOP04Fo.5. The second low-temperature 
phase, the 7'-polymorph with the short-range order of 
Na atoms (Sec. |III[ ), is more difficult to model due to 
the unresolved short-range order, although we may ex- 
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FIG. 18. (Color online) Magnetic layers in the crystal struc- 
ture of 7-Nai,5VOP04Fo.5 (left) and the respective spin lat- 
tice with six inequivalent couplings (right). The Ji — J2 pla- 
quette units are shaded. 



pect a similar scenario with FM couplings between near- 
est neighbors and AFM couplings between next-nearest 
neighbors. 

The LDA density of states (DOS) for 
7-Nai.5VOP04Fo.5 identifies oxygen and fluorine 
2p valence bands below —3 cV and vanadium 3d bands 
at the Fermi level (Fig. 19 1. The vanadium states 
show a crystal-field splitting into the t2g (below 0.7 eV) 
and eg (above 1.3 eV) sublevels, as expected for the 
octahedral coordination of a transition metal. Further 
on, the short bond to the axial oxygen atom (01) splits 
the t2g bands into the half-filled dxy states lying at 



the Fermi level and unoccupied dyz , dxz states lying 
between 0.2 eV and 0.7 eV (here, z aligns with the 
crystallographic c axis). The positions of the crystal-field 
levels justify our empirical conclusions on the symmetry 
of the magnetic orbital that induces leading couplings 
in the ab plane (Sec. IV A). The crystal-field splitting 

y resembles the positions of 
(Ref.[Sll), Ca(VO)2(P04)2 



in 7-Nai.5VOP04Fo.5 closel; 
the 3d sublevels in Pb2V309 
(Ref. 1 39 1), and other V"*"* oxides. 

The gapless (metallic) nature of the energy spectrum 
is a typical problem of LDA due to the underestimate 
of correlation effects in the V 3d shell. While the corre- 
lation effects can be introduced in a mean-field way via 
the LSDA-I-J7 method, the huge size of the unit cell (144 
atoms) tempts us to avoid this procedure and apply an 
LDA-based model approach instead. The latter is re- 
markably efficient for the evaluation of weak exchange 
couplings in compounds. 52j53, Using Wannier 

functions with proper orbital characters,!^ we fit 80 vana- 
dium 3d bands with a tight-binding model, and map the 
resulting hopping parameters (t) onto the multi-orbital 
Hubbard model with the effective on-site Coulomb repul- 
sion UeB and on-site Hund's coupling Jeff. The notably 
low t values (below 5 meV) compared to UcS — 4 eV 
identify the strongly correlated limit t <C [4ff so that the 
expression of the Kugel-Khomskii model for the exchange 




Energy (eV) 

FIG. 19. (Color online) LDA density of states 
Nai.5VOP04Fo.5. The Fermi level is at zero energy. 



for 7- 



couplings can be applied I^SISS] 

E 



(C/eff + A„)(C/cff + A„ - Jeff) ■ 



(7) 



Here, a denotes unoccupied d orbitals, txy are hoppings 
between the xy orbitals, txy^a are hoppings from the xy 
(half-filled) to a (empty) orbitals, and Aq are the crystal- 
field splittings between the xy and a orbitals. The first 
term is AFM superexchange, whereas the second term 
is F M and arises from the Hund's coupling i n the 3d 
shelLl^SIIIIlll Using U^^ = 4 eV and Jeff = 1 eVW^ we 
arrive at individual exchange couplings listed in Table|Vl] 

The weak exchange couplings in Nai.5VOP04Fo.5 are a 
severe challenge for computational methods. In contrast 
to PbZnVO(P04)2 and Pb2V309 (Refs. [15]and[S31), the 
calculated J's are exclusively AFM, despite the fits of the 
magnetic susceptibility (Sec. IV A[ ) reveal the presence of 
FM exchange. This problem is likely unavoidable and 
restricts the DFT results to qualitative conclusions that 
are nevertheless essential for understanding the system. 
The sizable AFM NNN couplings, compared to weaker 
NN couplings, support the FM Ji - AFM J2 scenario. 
The weakest NNN exchange J2 is about 70 % of J2 and 
J2 so that a spatial anisotropy of the diagonal couplings 
is moderate, as in Pb2VO(P04)2 and PbZnVO(P04)2 
(Refs. [T?l and By contrast, the difference between 
J" and J[ (Ji) is about 5 K, which is larger than the 
experimental | Ji| — 3.7 K. This suggests a sizable spatial 
anisotropy of the NN couplings. 

Owing to the complex structure of the spin lat- 
tice, the relation between the averaged couplings (Ji, 
J2) and individual exchanges is more intricate than in 
AA'VO(P04)2.'^ Since there are two Ji bonds yet sin- 
gle J[ and Jj" bonds per lattice site, one writes Ji = 
(2Jj + JI + J'i)/^- Applying this expression to the DFT 
estimates in Table [Vl| we arrive at J2 = 13.4 K that is 
twice larger than the experimental J2 ~ 6.6 K. A similar 
two-fold overestimate is found in AA'VO(P04)2 phos- 
phates: compare J2 = 18.4 K from DFT to J2 = 10.0 K 
from the experiment in PbZnVO(P04)2 (Ref. 25). This 
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TABLE VI. FM ( J^"^), AFM (J^^"^), and total (J) exchange 
integrals calculated with Eq. Q. The intralayer couplings are 
depicted in Fig. |18| whereas J±_ couples the neighboring layers 
(see Fig. [T]). 
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gives a strong support to the DFT results, and ascribes 
their apparent inaccuracy to systematic errors that do 
not alter qualitative trends. 

The notable reduction in J2 (6.5 K in Nai 5VOPO4F0.5 
vs. 9 — 10 K in PbZnVO(P04)2 and similar phosphates) 
can be traced back to longer bonds comprising the su- 
perexchange pathways. The relevant parameters are 
V-0 and 0-0 distances.f^Sl The latter are edges of the 
PO4 tetrahedra and stay nearly constant, whereas the 
former are flexible. In Nai^OP04Fo.5, the V-O dis- 
tances exceed 2.0 A (Table [lll| and contrast with shorter 
(1.95 — 2.0 A) distances for leading AFM couplings in 
AA'VO(P04)2 (Table IV in Ref.lll]). Nai.5VOP04Fo.5 
resembles Li2VOSi04 in terms of the diagonal couplings 
(J2 = 6.7 K and 6.3 K, respectively), despite the fact 
that the PO4 anion is smaller than Si04. 

While the nature of the diagonal couplings in 
Nai.5VOP04Fo.5 is readily comprehended, the evalua- 
tion of the NN couplings remains a challenge because 
of the computational inaccuracies of DFT. To make 
things even more complicated, the effects of the spa- 
tial anisotropy cannot be resolved by the available ex- 
perimental data on Nai.5VOP04Fo.5. The orthorhom- 
bic distortion of the FSlI^ favors a certain direction 
for the spin columns, stabilizes the columnar AFM state 
and, therefore, increases the saturation field H^. When 
the experimental Hg exceeds the prediction of the regu- 
lar FSL model, their difference approximates the magni- 
tude of the orthorhombic distortion for NN couplings.^- 
The spin lattice of Nai.5VOP04Fo.5, however, lacks any 
orthorhombic distortion and rather shows a complex 
tetragonal superstructure with three inequivalent cou- 
plings (Jl, J{, and J"). The couplings of each type run 
along both dimensions of the square lattice (Fig. 18 ) and 
do not stabilize the columnar AFM state, thereby no ef- 
fect on Hs should be expected. Presently, we are un- 
able to give an experimental-based estimate of the spatial 
anisotropy in Nai.5VOP04Fo.5. As long as the DFT re- 
sults are considered, appreciable anisotropy effects should 
be expected (Table VI ) . 

The last important feature of Nai.5VOP04Fo.5 is the 




20. (Color online) Vanadium dxy bands at the Fermi 
level (zero energy): LDA band structure (thin light lines) and 
the fit of the tight-binding model (thick dark lines). 



sizable interlayer coupling Jl. It is related to dimers 
of octahedra sharing the fluorine atom (right panel of 
Fig. [1]) , and couples each lattice site to one of the neigh- 
boring layers only. Since the magnetic saturation requires 
all spins to be parallel, the presence of the AFM inter- 
layer coupling increases the saturation field. Eq. ^ can 
be rewritten as 



/ioiJ. = (2Ji + 4J2 + J±)kB/igfiB) 



(8) 



to yield J± = 1.1 K, where we used Ji = —3.7 K, J2 — 
6.6 K, ^qHs = 15.4 T, and g = 1.95. This experimental 
estimate of J± compares well to the calculated value of 



1.7 K (Table VI) and provides additional justification for 
the reliability of our DFT results. 

The interlayer couphng in Nai.5VOP04Fo.5 amounts 
to 20 % of J2 and is larger than in any of the known 
FSL compounds (see Table VII I. Jj_ arises from the 
short V-V distance of 4.21 A, which is in fact the short- 
est V-V separation in the structure. There are compa- 
rable V-V separations of 4.45 - 4.55 A in Li2VOSi04 
and Li2VOGe04, but the lack of the bridging fluorine 
atom reduces the interlayer coupling to Jj^ ~ 0.2 K and 
drives these compounds to a perfectly 2D regime.'^^'^ 
The ^^'VO(P04)2 phosphates show even weaker inter- 
layer couphngs J_L ^ 0.1 kP^ owing to the V-V separa- 
tions of 8.5—10.0 A. The influence of the sizable interlayer 
coupling on the magnetism of Nai.5VOP04Fo.5 is further 
discussed in Sec. IVIBI 



VI. DISCUSSION 
A. Crystal structure 

The crystal structure of Nai.5VOP04Fo.5, as reported 
in Refs. [26] and [27l is deceptive in its simplicity. It dif- 
fers from our findings for Nai.5VOP04Fo.5 prepared by 
high-temperature solid-state method. The a-type struc- 
ture, previously found at RT, is observed above 500 K 
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FIG. 21. (Color online) The arrangement of Na atoms in different polymorphs of Nai.5VOP04Fo.5. Dark spheres denote filled 
Na positions, light spheres are partially filled Na positions, and open circles represent vacancies. 



only. The actual RT structure is marginally different, 
owing to the partial ordering of the Na atoms in the /3- 
polymorph. Low-temperature studies reveal the forma- 
tion of the complex superstructure in the 7-polymorph 
and the (3 ^ 7 + 7' phase separation. The apparent 
differences in the RT structure are probably related to a 
variety of preparation procedures that could lead to slight 
deviations in the sample composition, especially the 0/F 
ratio (see Sec. |l|. In our case, magnetization measure- 
ments combined with ESR and NMR local probes con- 
firm the oxidation state of -1-4 for the vanadium atoms 
and ensure the ordered arrangement of O and F. Since 
our samples are synthesized from stoichiometric mixtures 
of the reactants, the Nai.5VOP04Fo.5 composition is 
safely established. 

The course of the structural transformations 
in Nai.5VOP04Fo.5 is reminiscent of isotypic 
Nai.5AfP04Fi.5 fluorophosphates with = Fe, 

Cr, V, Al, and Ga.^^ These compounds accommodate 
trivalent M cations by substituting F for 01. This 
way, the M cation attains a nearly regular octahedral 
coordination with fluorines in the axial positions. 
High-temperature structures are always of the a-type, 
but low- temperature features are variable. When M — 
Fe and V^'^, the RT structure is of the /3-type, with 
a — v^asub, P4:2/mnm space group, and partially 
ordered Na atoms. The complete ordering of Na takes 
place in the orthorhombic 7-phase (a ~ 6 ~ 2asub) below 
RT. By contrast, the fluorophosphates with AI = Cr, Al, 
and Ga drop directly to a 7-type RT polymorph having 
a — 2asub, P4:2/mbc space group, and still disordered 
Na atoms. Below RT, these compounds presumably 
undergo another phase transition, but its nature has 
not been established. Nai.5VOP04Fo.5 complements 
this welter of polymorphs and transitions by the ordered 
tetragonal 7-phase and the /? -H> 7 + 7' phase separation 
below RT. 

Although the origin of the structural transforma- 
tions in Nai.5A/P04Fi 5 is far from being clear, elec- 
tronic effects should likely be ruled out. Indeed, mag- 



netic (Cr) and non-magnetic (Al, Ga) trivalent cations 
give rise to the same sequence of the transitions.!^ In 
Nai.5VOP04Fo.5, the effect is even more tangible. Weak 
exchange couplings manifest themselves below 10 — 20 K 
which is far below the flrst-order structural transition 
and phase separation around 250 K. In Fig. 21 we plot 



the arrangement of Na atoms in different polymorphs of 
Nai.5VOP04Fo.5. The partial filling of the two multiple- 
site positions in the a-polymorph is an indication of 
Na mobility. /3-Nai.5VOP04Fo.5 reveals clusters of or- 
dered (Nal) and partially disordered (Na2) atoms. In 
7-Nai.5VOP04Fo.5, we find a peculiar ordering pattern 
reminiscent of a ^-depleted square lattice. A similar 
pattern is observed in 7-Nai.5FeP04Fi.5, the only ex- 
ample of the complete cation ordering in this structure 
typcf^The structural transformations/phase separation 
are possibly caused by electrostatic interactions between 
Na and F atoms together with an apparent difficulty of 
finding a symmetric configuration for the |-filled "square 
lattice" of Na. 



A notable feature of the phase separation in 
Nai.5VOP04Fo.5 is the lack of any visible compositional 
changes. When electronic degrees of freedom are inac- 
tive, phase separation is typically driven by a size mis- 
match and results in different compositions of the two 
phases.!^ In Nai.5VOP04Fo.5, the fully ordered struc- 
ture of the 7-phase implies the ideal Nai.5VOP04Fo.5 
composition and imposes the same composition for the 
a-phase, as further confirmed by the unique oxidation 
state of -1-4 for vanadium. Neither ESR nor NMR re- 
veal any differences between the two phases coexisting at 
low temperatures. Another surprising observation is the 
macroscopic scale of the separated phases: both 7- and 
7'-polymorphs show narrow XRD reflections correspond- 
ing to well-deflned crystallites rather than nanoscale do- 
mains. These unusual features suggest that a further 
study of the phase separation in Nai 5VOPO4F0.5 could 
be insightful. 
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B. Magnetic properties 

Our detailed study of the magnetic properties 
and electronic structure confirms the assignment of 
Nai.5VOP04Fo.5 to the quasi-2D FSL-type spin model. 
The ultimate goal of finding an ideal FSL with FM Ji 
is, however, defeated by the complexity of structural 
transformations and low-temperature phase separation. 
The structural changes are not caused by the magnetic 
frustration (Sec. |VI"A' ), while the magnetism, in turn, is 
marginally influenced by the spin lattice distortion. Sim- 
ilar to AA'VO(P04)2, thermodynamic properties con- 
form to the regular FSL model. A closer and micro- 
scopic look, based on DFT, detects sizable deviations 
from the regular model, although experimental data give 
little information on such effects. In contrast to our ear- 
lier conjecture in Ref. \Wl the underestimate of the satu- 
ration field should be ascribed to the interlayer coupling 
(Sec.|v]). 

Despite seeming to be overcomplicated at first glance, 
the spin lattice comprising six inequivalent intralayer 



couplings (right panel of Fig. 18 ) expands the family of 
FSL-type models realized in inorganic compounds. The 
AA'VO(P04)2 phosphates reveal a somewhat more sim- 
ple, orthorhombic distortion of the FSL.^'* By contrast, 
the tetragonal superstructure in 7-Nai.5VOP04Fo.5 
splits the square lattice into plaquettes (Fig. 18) and 



therefore relates to a different branch of theoretical re- 
sults. The models of tetramerized square lattices, also 
known as square lattices with plaquette structure, ex- 
hibit a quantum phase transition between the long-range- 
ordered and gapped ground states but so far lack any pro- 
totype materials. A cation substitution/intercalation in 
Nai.5VOP04Fo.5 could be a way to obtain further com- 
pounds of this type, although the possible problem of the 
phase separation should not be overlooked. 

Apart from the spatial anisotropy, Nai.5VOP04Fo.5 
features an unusually strong interlayer coupling. In Ta- 
ble VII we list some basic parameters of various FSL 



compounds with columnar AFM ground state: experi- 
mental values of Ji and J2 , thermodynamic energy scale 
Jc ~ i/Jf -I- J|, DFT estimates of J±, and experimen- 
tal temperatures of the magnetic ordering (T/v). To ac- 
count for the single interlayer coupling per lattice site in 
Nai.5VOP04Fo.5, we take one half of the experimental 
J±. The comparison identifies the effects of the frustra- 
tion and interlayer couplings on T/v- Based on Tjv /Jc, the 
compounds can be divided into three groups: strongly 
frustrated {a ~ —1, T^/Jc — 0.2), moderately frus- 
trated [a ~ —2, Tisi/Jc — 0.35), and weakly frustrated 
[a > 5, Tjy/Jc — 0.45 — 0.50) systems. An order of mag- 
nitude difference in J±/Jc between Nai.5VOP04Fo.5 and 
PbZnVO(P04)2 has no appreciable effect on T/v. 

Despite the weak influence on T/v, J± likely affects 
other features related to the magnetic ordering. In 
Fig. |22[ we compare specific-heat transition anomalies for 
Li2VOSi04, Nai.5VOP04Fo.5, and Pb2VO(P04)2- The 
area under the Cp/T curve is a measure of the transi- 
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FIG. 22. (Color online) Zero-field specific heat at the mag- 
netic ordering transition in Nai.5VOP04Fo.5 (powder sam- 
ple), Li2VOSi04 (powder sample), and Pb2VO(P04)2 (single 
crystal).''^ The data for the last two compounds are from 
Ref.lISl 



tion entropy Pb2VO(P04)2 shows the smaUest 5, 
below 0.01 J mol~^ K~^. A similar frustration ratio of 
a ~ —2 leads to a much larger S ~ 0.065 J mol^^ K^^ 
in Nai.5VOP04Fo.5, while Li2VOSi04 reveals an even 
broader anomaly with S' ~ 0.11 J mol~^ K~^. Although 
numerical estimates of S are tentative, the qualitative 
trend is robust and identifies the strong dependence of 
5 on Jj^. The reduction in Jj_ shifts the entropy to the 
broad maximum above Tjv and reduces the transition 
entropy, which determines the magnitude of the tran- 
sition anomaly. At weak J_\_, the anomaly ultimately 
transforms into a kink, especially for powder samples, 
see the data in Refs. [13 [Ml and [251 A similar trend for 
the transition entropy depending on Jl has been pro- 
posed in a theoretical studjH^ of a simplified spin model 
with non- frustrated square lattices coupled by Jj_. We 
mention that the 3D derivative of the FSL model has 
been addressed in recent theoretical studies focused on 
ground-state properties.!^ Our results call for a further 
investigation of finite-temperature behavior: in particu- 
lar, the evaluation of Neel temperature and temperature 
dependence of the specific heat. 

Finally, Nai.5VOP04Fo.5 can be a starting point to ex- 
plore FSL-type systems with different values of the spin. 
Experimental work has been mostly restricted to the case 
of spin-i,!^ although theoretical results propose pecu- 
liar features of systems with larger spin.*^ The family of 
Nai.5MP04Fi.5 fluorophosphated^ should retain the ba- 
sic FSL-type geometry with FM couplings between near- 
est neighbors and AFM couplings between next-nearest 
neighbors in the ah plane. Trivalent M cations represent 
spin-1 (V+^), spin-| (Cr+^), and even spin-| (Fe+^) so 
that the effect of the spin value on the FSL model could 
be explored experimentally. 

In summary, we have presented a detailed study of 
the crystal structure, magnetic behavior, and electronic 
structure of Nai 5VOPO4F0.5, a spin-^ FSL-type com- 
pound. We observed essential structural changes, includ- 
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TABLE VII. Comparison of FSL compounds with the columnar AFM ground state: exchange 
couplings Ji and J2 (in K); thermodynamic energy scale Jc ~ \/ Jj + J| (in K); the frustration 
ratio Q = J2/J1; the interlayer coupling Jx; and the magnetic ordering temperature Tjv. 



Compound 


Ji 


J2 


a 


Jc 


J±/Jc 


Tjv / Jc 


Ref. 




— 0.0 




n Q 






U. Z i 




SrZnVO(P04)2 


-8.3 


8.9 


-1.1 


12.2 


< 0.01 


0.22 




BaZnVO(P04)2 


-5.0 


9.3 


-1.9 


10.6 


< 0.01 


0.36 


\m 


Pb2VO(P04)2 


-5.2 
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-1.8 


10.7 


< 0.01 


0.33 
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PbZnVO(P04)2 


-5.2 


10.0 


-1.9 


11.3 


< 0.01 


0.35 
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Nai.5VOP04Fo.5 


-3.7 


6.6 


-1.8 


7.6 


0.07 


0.34 


This work 


Li2VOSi04 


0.6 


6.3 


+10.5 


6.3 


0.03 


0.44 




Li2VOGe04 


0.8 


4.1 


+5.1 


4.2 


0.05 


0.50 





ing the low-temperature phase separation, that prevent 
Nai.5VOP04Fo.5 from being a simple model system. A 
complete experimental evaluation of the spin model for 
this compound is a formidable challenge. Despite this 
fact, Nai.5VOP04Fo.5 is the only experimental example 
of the tetramerized (plaquette-structure) square lattice 
and bears a relation to a group of spin models which so far 
lacked any prototype materials. Our experimental data 
conform to the regular FSL model supplied with a siz- 
able interlayer coupling J± of about 10 % of the effective 
intralayer exchange. This interlayer coupling marginally 
affects the magnetic ordering temperature, yet the transi- 
tion anomaly in the specific heat is notably increased. In 
general, Nai.5VOP04Fo.5 and isostructural compounds 
establish several promising connections to FSL-type spin 
models. This should outweigh the apparent structural 
complexity of these materials and stimulate further ex- 



perimental investigation. 
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FIG. SI. Structure refinement of a-Nai.5VOP04Fo.5 at 560 K: experimental (symbols), 
calculated (green line), and difference (black line) patterns. Ticks show reflection positions. 
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FIG. S2. Structure refinement of /3-Nai.5VOP04Fo.5 at RT: experimental (symbols), calcu- 
lated (green line), and difference (black line) patterns. Ticks show reflection positions. 
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FIG. S3. Structure refinement at 150 K: experimental (symbols), calculated (green line), and 
difference (black line) patterns. Upper and lower sets of ticks denote reflections of the 7'- 
and 7-polymorphs, respectively. 
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FIG. S4. Specific heat of Nai.5VOP04Fo.5 measured in zero magnetic field. 



